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SOLAR  AMDSPHERIC  DYNAMICS 

I.  Introduction 

The  goal  of  this  research  Is  to  study  the  drl\riJTg  of  the  solar  atmospheric 
oscillations  and  their  role  In  heating  the  chromosphere  and  corona.  During 
the  first  year  of  this  contract  the  wave  vector  surfaces  of  I'fegneto-Acoustic- 
Gravity  (M.A.G.)  waves  were  studied  and  we  derived  modal  equations  for  conpress- 
Ible  three-dimensional  fluid  flow  that  will  enable  the  propagation  and  dissi- 
pation of  i^.A.G.  waves  to  be  studied.  During  the  second  year  a ccmputer  code 
to  solve  these  modal  equations  was  written  and  testing  on  it  begun.  In  addi- 
tion, data  from  Skylab  was  analyzed  for  acoustic  pulses  in  the  transition 
region  and  the  steady  flow  of  a solar  wind  through  a magnetic  flux  tube  diverg- 
ing more  rapidly  than  radially  was  calculated.  The  contract  supported  Dr. 

S.H.  Hsieh  as  a research  fellow  and  Mr.  Robert  Wolff  as  a graduate  student 
working  on  the  solution  of  the  modal  equations.  So  far  under  this  contract 
one  paper  has  been  published  on  Magneto-Gravity  Waves  (Schwarts  and  Stein, 

1975),  and  two  others  are  now  being  prepared.  One  on  "Magnetc-Acoustic-Gravity 
Waves"  and  the  other  on  "Flow  Through  a Diverging  i-lagnetic  Flux  Tube". 

II.  Wave  Vector  Surfaces  of  .^^agneto-Acoustic-Gravity  Waves 

I have  studied  the  wave  vector  surfaces  of  Kaigneto-Acoustic -Gravity  (MAG) 
waves  in  the  WKE  limit,  that  is  using  the  local  dispersion  relation.  The  -/iarje 
vector  surface  is  the  surface  k(u))  in  k-space  on  which  the  local  dispersion 
relation  is  satisfied  for  fixed  frequency  w. 
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There  are  three  wave  modes.  One  is  the  usual  Alfven  mode.  It's  wave 
vector  surface  consists  of  two  planes  perpendicular  to  the  magnetic  field 
and  a distance  o)/acos0  from  the  origin.  We  call  these  the  Alfven  planes. 

At  lilgh  frequencies,  w>N  , gravity  is  unimportant.  Here  N =rg/2s  is 
the  acoustic  cutoff  frequency,  where  g is  the  acceleration  of  gravity,  s is 
the  sound  speed  and  r the  ratio  of  specific  heats  at  constant  pressure  and 
volume.  In  this  case,  the  MAG  waves  betiave  as  Magneto-hydrodynamic  (MHD) 
waves.  In  weak  magnetic  fields  the  fast  mode  is  an  acoustic  mode  and  the 
slow  mode  a magnetic  mode,  while  in  strong  magnetic  fields  the  fast  mode  is 
a magnetic  pressure  mode  and  the  slow  mode  an  acoustic  mode  propagating 
along  the  magnetic  field. 

At  low  frequencies,  one  mode  has  a wave  vector  surface  that  is  approxi- 
mately a plane  perpendicular  to  the  magnetic  field.  This  wave  vector  surface 
lies  everywhere  closer  to  the  origin  than  the  Alfven  plane  and  has  an  outward 
bunp  that  Just  touches  the  Alfven  plane  for  vertically  propagating  waves 
(kj^orlz”^^'  the  wave  vector  surface  approaches  a plane  that 

passes  through  the  origin. 

The  other  mode's  wave  vector  surface  is  similar  to  the  acoustic-gravity 
'v^^ave  hyperboloid  close  to  the  origin  and  also  becomes  a plane,  perpendicular 
to  the  magnetic  field  far  away  from  the  origin.  As  the  magnetic  field  strength 
increases,  the  hyperboloid  gets  squashed  between  the  Alfven  planes  and  dis- 
appears. In  the  regions  of  large  neither  wave  is  a normal  propagating 

wave.  For  both  their  potential  energy  becomes  far  greater  than  their  kinetic 
energy.  In  the  weak  field  case  the  magnetic  potential  energy  is  small  for 
both  modes  when  large,  while  in  the  strong  fiexd  case  one  mode  has 

predominantly  magnetic  potential  energy  and  the  other  has  none. 


Part  of  this  work  was  published  by  Schwartz  and  Stein  (1975)  and  part 
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appeared  in  a thesis  of  Nelson  Hartunian  (1975)  from  Brandeis  University. 
The  remainder  is  now  being  written  up  for  publication.  A clearer  picture 
of  MAG  waves  will  however  require  the  solution  of  the  initial  value  problem 
either  with  the  modal  or  a truly  3-dime.nslonal  code. 


III.  fiodal  Eauations  for  Compressible  Fluid  Flows 


Non-linear  calculations  are  simple  in  one  spatial  dimension  (Richtmyer 
and  rtorton,  1969),  but  conpllcated  in  two  or  three.  However,  gravity  waves 
are  inherently  two-dimensional  and  magneto-acoustic-gravity  waves  are  inher- 
ently three-dimensional.  The  modal  approach  reduces  the  three-dimensional 


problem  to  a one-dimensional  problem  by  imposir.g  a horizontal  structure  on 
the  motion.  We  separate  all  the  fluid  variables  into  mean  and  fluctuating 
parts.  The  mean  part  is  urMform  on  horizontal  planes  and  the  fluctuating 
part  has  a periodic  hexagonal  horizontal  structure  of  some  given  length  scale. 
We  then  have  a system  of  Eulerian  partial  differential  equations  for  the  mean 
and  fluctuatlr.g  variables.  These  equations  were  derived  in  collaboration 
’/d.th  H.S.  Hsieh  and  Robert  Wolff. 

To  solve  these  equations  we  write  them  in  findte  difference  "conservation" 


form  and  use  the  Richtmyer  two-step  Lax-Wendroff  method  (Plichtmyer  and  Thorton, 


j 
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The  Lax-Wendroff  scheme  is  dissipative  and  handles  shocks  well,  but  it  also 
dairps  waves  with  wavelengths  less  than  20  zones  moderately.  It  is  also  subject 
to  some  non-linear  instabilities.  To  suppress  these  instabilities  we  introduce 
additional  dissipation  with  the  "flux  corrected  transport"  scheme  developed  by 
J.  Boris  (Book,  Boris,  Hain  1975).  The  Lax-Wendroff  scheme  was  tested  for  the 
one-dimensional  full  fluid  equations  on  the  shock  tube  problem  and  a small  am- 
plitude piston  driven  wave  in  a uniform  medium.  Tiie  shock  tube  problem  is  an 
especially  sensitive  test  because  it  involves  three  types  of  waves  (shock,  contact 
discontinuity  and  simple  wave)  l^rlth  known  speeds  and  amplitudes,  and  is  independ- 
dent  of  boundary  conditions. 

The  boundary  conditions  are  more  difficult  for  the  modal  equations  than 
for  Lagrangian  one-dimjsnsional  equations.  The  Rlchtmyer  Lagrangian  schem.e  is  a 
leapfrog  scheme  with  the  variables  staggered  in  space  and  time.  Hence  the 
boundary  conditions  need  only  specify  the  velocity.  The  Lax-Wendroff  scheme 
defines  all  variables  at  the  same  space  and  time.  Hence  we  must  specify  all 
the  variables  at  the  boundary.  To  do  this  we  need  to  know  the  solution  at  one 
point,  but  in  general  we  do  not.  In  our  case,  we  can  use  polarization  relations 
for  small  amplitude  waves  to  specify  all  variables  at  the  inflow  boundary.  At 
the  outflow  boundary,  we  can  use  essentially  the  same  method  as  in  the  Lagran- 
gian scheme.  For  any  variable  f at  the  boundary  Xj, 


These  boundary  conditions  have  been  coded  for  acoustic  waves  in  a uniform  medi- 
um and  the  outflow  boundary  conditions  for  a stratified  atmosphere  as  well. 

It  remai.ns  to  code  the  Lnflow  boundary  conditions  for  acoustic-gravity  -fi3.ves. 

J 
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A conputer  code  for  the  modal  equations  (without  magnetic  field  or 
radiation)  is  currently  being  tested  in  collaboration  with  Robert  Wolff.  The 
testing  program  is  as  follows: 

A.  To  test  non-linear  coupling,  consider  plane-parallel  (vertical) 
motion  in  a uniform  medium  for 

(1)  shock  tube 

(2)  piston  driven  sinusoidal  waves  Including  shock  formation 

B.  To  test  horizontal  variation,  consider  a piston  driven  wave  in  a uni- 
form medium.  This  corresponds  to  waves  propagating  at  an  angle  to 
the  vertical  direction.  Test: 

(1)  increase  in  vertical  ^'ra.velength  with  decreasing  horizontal 
scale 

(2)  velocity  is  longitudinal 

(3)  reduced  height  of  shock  formation 

(4)  boundary  condition  for  non-vertically  propagating  wave 

C.  To  test  vertical  motion  in  a stratified  atmosphere,  consider  piston 
driven  vertical  acoustic  waves  v\rith  frequency  above  and  below  the 
critical  frequency  and  compare  with  the  one-dimensional  Lagrangian 
calculations  usirig  the  full  fluid  equations . 

Once  the  code  has  been  fully  tested  for  acoustic -gravity  vjaves,  conduction 
and  radiative  losses  will  be  added  and  a heating  model  of  the  solar  chromosphere 
and  corona  constructed.  Afterwards  magnetic  fields  and  grey  radiative  transfer 
will  be  added.  Simultaneously,  grey  radiative  transfer  vrtll  be  added  to  the 
one-dimensional  Lagrangian  code.  The  modified  code  \vill  be  used  to  construct  a 
dynamic  model  solar  atmosphere.  This  will  be  used  as  an  initial  model  for  a 
non-equilibrium,  non-grey  time  dependent  model,  which  vri.ll  be  used  to  analyze 
changes  observed  in  line  profile  due  to  waves  propagating  through  the  atmosphere 
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in  collaboration  with  Richard  Klein  (Kitt  Peak  National  Observatory)  and 
Richard  Shine  (University  of  Colorado).  Richard  Shine  is  currently  calculating 
Ha  emission  from  the  accretion  shock  found  in  the  calculations  by  Klein,  Stein, 
and  Kalkofen  (1977)  of  a radiating  shock  propagating  through  an  A star  atmos- 
phere . 

The  Lax-Wendroff  scheme  is  also  being  used  in  the  development  of  a tv;o- 
and  three-dimensional  code  in  collaboration  with  Richard  Klein  and  Lawrence 
Auer  (High  Altitude  Observatory,  National  Center  for  Atmospheric  Research). 

Such  a code  will  be  able  to  solve  truly  three-dimensional  magneto-hydrodyramic 
problems.  We  have  so  far  made  tests  on  two  two-dimensional  problems  — a 
symmetric  shock  and  gravitational  accretion  flow  into  a density  perturbation. 

We  have  tested  thnee  schemes:  MacCormack  (1971),  Lax-Wendroff  and  Durner  Cell 

(Book,  Boris,  and  Hain,  1975).  So  far  only  the  MacCormack  scheme  is  conpatible 
ivith  our  Poisson  equation  solver  for  the  gravitational  potential,  but  all  tliree 
work  for  the  shock  problem.  We  are  now  working  on  the  data  manage.ment  problem 
for  three-dimensional  calculations  and  on  a faster  Poisson  solver. 


IV.  Wave  Pulses  in  the  Transition  Region 


Vertically  propagating  acoustic  pulses  develop  a secondary  crest  in  their 
'flake  (Figure  1).  The  second  peak  arises  from  the  resonari  wake  behind  pulses 
in  a stratified  atmosphere  (Lamb,  19^5;  Schmidt  and  Zirker,  1963;  Kato,  1966; 
Stein  and  Schwartz,  1972).  The  separation  betv/een  the  peaks  is  close  to  the 
acoustic  cutoff  period  of  the  atmosphere.  This  two  peaked  wave  form  may  be 
used  as  a signature  for  acoustic  pulses.  Obser'/atlons  from  Sky lab  show  inten- 
sity fluctuations  in  transition  region  lines  that  occasionally  displayed  such 
a double  peaked  structure  (Figure  2)  (Vemasza  et  al  1975).  Vemazza  and  I 
have  made  a systematic  search  for  such  double  pealced  pulses,  and  performed  some 


; calculations  on  pulse  propagation  to  determine  the  relation  between  density 

and  velocity  fluctuations. 

If  acoustic  pulses  with  their  double  peaked  structure  are  present,  they 
' should  become  more  prominant  as  we  look  at  stronger  pulses.  Weak  pulses  will 

be  common  and  precede  other  weak  pulses  by  various  time  intervals  without  any 
. causal  connection.  However,  strong  pulses  are  rare,  hence  two  such  pulses 

following  one  another  closely  in  time  are  likely  to  be  causally  connected.  To 
search  for  the  presence  of  acoustic  pulses  in  the  Skylab  data,  we  made  histo- 
grams of  the  number  of  pulses  seen  in  CIII  and  OIV  intensity  greater  than  a 
' given  strength  versus  the  interval  between  succeeding  pulses  (Figures  3 and  4). 

'when  all  pulses  are  included,  the  number  of  pulses  at  any  given  separation 
decreases  nearly  monotonlcally  with  increasir.g  pulse  separation  from  a maxi- 
mum at  the  minimum  resolvable  separation.  As  the  pulses  included  are  restric- 
ted to  those  of  greater  and  greater  strength  a decided  maximum  develops  in  the 
range  of  100-175  seconds  separation.  This  behavior  is  clearly  seen  i.n  the  Or/ 
data.  The  data  for  CIII,  which  is  formed  at  9x10**  ~K,  is  more  complicated.  In 
the  case  of  CIII,  even  for  the  weakest  pulses,  there  is  a maximum  near  150  sec. 
separation,  and  this  maximum  persists  as  the  pulse  strength  increases. 

A separation  of  150  sec.  corresponds  to  the  acoustic  cutoff  period  In  an 
isothermal  atmosphere  of  3000  °K.  Note  there  is  also  a secondary  maximum  in 
i both  the  CIII  and  Or/  data  at  a separation  of  250  sec . , which  corresponds  to 

the  acoustic  cutoff  period  in  an  isothermal  atmosphere  of  8000  °K. 

Next  we  looked  at  why  ivaves  are  difficult  to  detect  in  intensity  fluctuations 
in  the  transition  region.  We  were  especially  interested  in  why  the  300s  oscilla- 
tions stop  beiTig  seen.  As  a simple  preliminary  calculation  we  kept  the  tempera- 
ture of  each  fluid  blob  fixed  at  its  initial  value  as  the  blob  moved  in  response 
to  waves.  This  case  corresponds  to  infinitely  rapid  radiative  energy  adjust- 
ment flow.  The  results  are  shown  in  Figure  5.  We  see  that  when  an  acoustic 

L-i 
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wave  enters  the  transition  region  its  velocity  ainplitude  continues  to  Increase 
with  height,  but  its  density  fluctuation  amplitude  Sp/p^  decreases  by  a factor 
of  3-  Hence  the  intensity  fluctuation,  which  is  proportional  to  density  squared, 
will  decrease  by  an  order  of  magnitude. 


V.  Flux  Tube  Flow 

Kopp  and  Holzer  (1976)  investigated  the  expansion  of  the  solar  'wind  along 
a tube  thiat  widens  out  faster  than  radially.  Such  a widening  tube  may  occur 
in  the  expansion  of  the  solar  wind  along  open  magnetic  field  lines.  They 
found  that  when  the  tube,  along  which  the  fl'uld  flows,  widens  by  more  than 
a critical  amount  the  flow  develops  two  sonic  points.  In  the  standard  case 
of  a radially  diverging  flow,  the  fluid  accelerates  and  becomes  supersonic  at 
a sonic  point  located  at  several  solar  radii.  When  the  flow  diverges,  more  rap- 
idly than  radially  close  to  the  sun,  a second  sonic  point  develops  in  the  region 
of  rapid  divergence.  Kopp  and  Holzer  investigated  the  time-dependent  transient 
behavior  of  the  solar  wind  when  the  flux  tube  along  which  it  is  flowing  suddenly 
'.'d.dens  by  more  than  the  critical  amount.  They  found  that  the  wind  accelerates 
and  becomes  supersonic  at  the  inner  critical  point,  shocks  back  to  the  subsonic 
solution,  then  accelerates  and  becomes  supersorJ-c  again  at  the  outer  critical 
point.  The  shock  is  facing  in'ward,  but  is  dragged  outivard  by  the  expanding 
wind.  In  the  course  of  several  hours  the  shock  is  dragged  cut  from  the  Lnner 
critical  point  to  the  outer  critical  point.  Thereafter,  the  flow  becomes  super- 
sonic at  the  inner  critical  point  and  remains  supersonic  throughout,  rejoining 
the  standard  solution  beyond  the  outer  critical  point. 

I repeated  these  calculations  using  an  Initial  value  Lagrangian  finite 
difference  code,  instead  of  the  Eulerlan  code  used  by  Kopp  and  Holzer.  I con- 
sidered the  same  case  as  they  had,  an  isothermal  atmosphere  of  pure  .hj.'drogen  at 
a tenperature  of  1.3  x 10°  °K.  The  area  of  the  flux  tube  'was  widened  from  a 
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radial  tube  to  a tube  which  diverged  by  a factor  4,5  greater  than  radial  at 
half  a solar  radius  above  the  surface.  The  formula  for  the  diverging  flux 
tube  was  the  same  as  used  by  Kopp  and  Holzer: 

A(r)/Ag  = (r/Rg)^  [f^  + f +1] 

where 

f = 1 - (f  -l)e''^o~^i^'^'^, 

1 max  ’ 

Contrary  to  their  results,  I found  the  shock  that  developed  outside  the  inner 
critical  point  remained  stationary  and  was  not  dragged  outward  by  the  wind 
(Figure  6) , 

Inside  the  shock  the  wind  speed  was  u^  = 231  km/s  = 1,33)  and  outside 
it  was  u^  = 104  km/s  (M^  = 0.6),  The  sound  speed  was  s^^  = 173  km/s.  The  vel- 
ocity Jump  across  the  shock 


was  found  from 

2 

y"  - [2  + (-5- 5-)  ] y + 1 = 0 

to  be  0.487.  The  shock  strer,gth  ’^as 


1/y  = 2.05, 

the  shock  velocity  ^/ra.s 


= - 17  km/s. 
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and  the  shock  Mach  number  was 
Mg  = y“^  = 1.43. 

Thus  the  shock  propagated  inward  very  slowly.  This  is  a stable  steady  state 
situation,  since  as  the  shock  moves  closer  to  the  critical  point  its  strength 
(1/y)  decreases  and  it  gets  dragged  back  out  to  where  it  was. 

The  Richtmyer  type  Lagrangian  code  I used  had  an  artificial  viscosity  which 
spread  the  shock  over  two  zones,  and  weakened  it  slightly.  The  real  shock  would 
therefore  be  stable  as  well,  and  would  be  slightly  further  in,  closer  to  the 
inner  critical  point.  The  reason  for  the  discrepancy  between  our  results  and 
those  of  Kopp  and  Holzer  is  not  known.  It  may  be  that  their  Eulerian  code 
dairped  their  shock  more  heavily  tbian  my  artificial  viscosity,  so  that  their 
shock  was  weaker,  and  hence  got  dragged  outward  by  the  wind. 

These  results  are  presently  being  written  up  as  a brief  research  note. 
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Figure  Captions 

Relative  intensity  fluct'uations  as  a function  of  tine.  The  relative 
intensity  is  shown  for  six  ILnes  formed  in  the  region  from  the  bottom 
of  the  transition  zone  (Ly  “)  to  the  low  corona  (Mg  x).  Intensity  is 
normalized  so  maximum  intensity  is  one  for  each  line  for  each  time 
interval.  Tenperatures  at  the  right  are  the  approximate  formation 
taiperature  of  the  lines. 

Velocity  as  a function  of  time  for  an  acoustic  '^fave  pulse  traveling 

vertically  upward  through  a model  solar  atmosphere.  The  velocity 

1/4 

amplitude  is  scaled  by  (density) 

Number  of  pulses  observed  in  0 IV  (x  554)  as  a function  of  the  time 
delay  between  successive  pulses.  Data  from  Harvard  UV  experiment 
on  sky lab.  Each  curve  is  for  pulses  greater  than  indicated  strength, 
as  measured  by  ratio  of  maximum  to  base  intensity. 

Similar  to  Figure  3,  but  for  C III  (x  977). 

Velocity  and  density  amplitude  of  pulse  at  bottom  (T=2xlO‘*°K)  and 
middle  (T=2xlO^°K)  of  transition  region.  Note  decrease  in  density 
but  not  velocity  amplitude  of  wave  in  transition  region. 

Mach  number  and  density  vs.  height  for  steady  flow  through  a ver- 
tical magnetic  flux  tube  that  diverges  faster  than  radially.  The 
flux  tube  spreads  by  a factor  4.5  faster  than  radially  at  heights 


near  Rq/2  above  the  photosphere. 
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